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ABSTRACT:  The  superionic  conductor  Cu2_,;Se  has  been  shown 
to  be  a  promising  thermoelectric  at  higher  temperatures  because  of 
very  low  lattice  thermal  conductivities,  attributed  to  the  liquid-like 
mobility  of  copper  ions  in  the  superionic  phase.  In  this  work,  we 
present  the  potential  of  copper  selenide  to  achieve  a  high  figure  of 
merit  at  room  temperature,  if  the  intrinsically  high  hole  carrier 
concentration  can  be  reduced.  Using  bromine  as  a  dopant,  we  show 
that  reducing  the  charge  carrier  concentration  in  Cu2_^Se  is  in  fact 
possible.  Furthermore,  we  provide  profound  insight  into  the 
complex  defect  chemistry  of  bromine  doped  Cu2_aSe  via  various 
analytical  methods  and  investigate  the  consequential  influences  on 
the  thermoelectric  transport  properties.  Here,  we  show,  for  the  first 
time,  the  effect  of  copper  vacancy  formation  as  compensating 
defects  when  moving  the  Fermi  level  closer  to  the  valence  band 

edge.  These  compensating  defects  provide  an  explanation  for  the  often  seen  doping  inefficiencies  in  thermoelectrics  via  defect 
chemistry  and  guide  further  progress  in  the  development  of  new  thermoelectric  materials. 


■  INTRODUCTION 

On  the  basis  of  Seebeck  and  Peltier  effects,  thermoelectric 
materials  either  convert  thermal  to  electrical  energy  or  vice 
versa,  providing  a  perspective  for  power  generation  and 
refrigeration  applications,  respectively.  The  thermoelectric 
efficiency  is  governed  by  the  thermoelectric  figure  of  merit 
zT  =  S1(pk)~1T,1'~  and  the  past  decade  has  seen  the  discovery 
of  many  suitable  materials  for  high  temperature  power 
generation.  One  of  these  materials  is  Cu2_^Se,  which  exhibits 
a  maximum  zT  of  1.6  at  1000  K,  making  it  a  competitive  and 
cheap  alternative  for  PbTe,  the  leading  thermoelectric  material 
in  the  hundreds  of  kelvin  temperature  range.  ’  Because  of  a 
small,  intrinsic  copper  deficiency  <5,  copper  selenide  (Cu2_^Se) 
is  a  p-type  semiconductor  with  a  band  gap  of  1.23  eV.+  Above 
410  K,  copper  selenide  transforms  to  a  Cu+  conducting  phase, 
wherein  the  copper  ions  are  mobile  and  described  as  “liquid 
like”  because  their  diffusion  coefficient  is  about  10-5  cm2s_1, 
comparable  to  the  value  for  water  molecules  in  liquid  water.’’0 
This  superionic  state  results  in  extremely  low  lattice  thermal 
conductivity  values  as  low  as  0.4  W(Rm)-1  in  the  superionic 
phase  and  in  turn  a  high  figure  of  merit.  Starting  with  the 
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superionic  Cu2_^Se,  many  other  superionic  copper  and  silver 
selenides  have  been  found  to  be  interesting  thermoelectrics;  for 
instance,  the  argyrodites  Cu7PSe6,7  Ag2Se,0-10  CuAgSe,11’1" 
and  quaternary  copper  chalcogenides  C^MM'Qt  (M  =  Zn, 
Fe,...;  M'  =  Zn,  Ge,  Q_=  S,  Se,  Te)  which  also  exhibit  interstitial 
copper  ions  and  even  show  band  convergence.1’  1  While 
copper  selenide  has  potential  as  a  thermoelectric  material  at 
high  temperatures,  its  ionic  conductivity  has  caused  long-term 
stability  problems  for  use  in  thermoelectric  generators.  °  This 
spurs  interest  in  the  properties  of  copper  selenide  below  the 
superionic  phase  transition  at  410  K,  where  the  ion  migration 
would  not  be  of  concern.  In  this  work,  we  investigate  the 
potential  for  copper  selenide  to  achieve  high  zT  values  at  room 
temperature.  Using  bromine  as  a  dopant,  we  show  the  effect  of 
carrier  scattering  on  the  thermoelectric  quality  factor  and  how 
the  dopant  bromine  itself  prevents  possibly  high  figures  of 
merit  at  room  temperature.  Furthermore,  we  explore  the 
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underlying  principle  for  doping  inefficiencies  in  copper 
selenides  via  the  intrinsic  defect  chemistry  of  this  material. 
We  show  that  doping  not  only  introduces  carriers  into  a 
semiconductor  but  also  changes  the  overall  stoichiometry  of  the 
host  material.  In  the  case  of  Cu2_aSe,  substitution  of  Se  with  Br 
introduces  compensating  defects  in  the  form  of  additional  Cu 
vacancies,  which  effectively  pin  the  Fermi  level  and  lead  to  a 
lower  doping  efficiency.  In  the  following  sections,  the  notation 
Cu2.ySe1_vBrl.  with  y  >  8  will  be  used  to  take  account  of  the 
additional  defects  in  the  Br  doped  samples. 

■  EXPERIMENTAL  SECTION 

Synthesis.  Bulk  samples  of  polycrystalline  Cu2.ySe1_;tBrl  with  x  =  0, 
0.02,  0.03,  0.04,  and  0.05  were  prepared  by  melting  and  annealing 
techniques  using  powders  of  Cu  (Alfa  Aesar,  99.999%)  and  Se  (Alfa 
Aesar,  99.999%),  as  well  as  CuBr  (Sigma,  99.999%).  In  order  to  ensure 
the  absence  of  copper  oxides  in  the  reaction  mixture,  the  copper 
powder  was  reduced  for  several  hours  at  525  K  under  N2/H2  (95:5) 
flow  in  a  tube  furnace  using  a  corundum  boat  as  the  reaction  vessel. 
Phase  purity  of  the  starting  materials  was  verified  via  X-ray  diffraction, 
and  all  synthetic  procedures  were  carried  out  in  a  N2  drybox.  The 
synthesis  was  performed  in  evacuated  quartz  ampules,  which  were 
dried  at  1073  K  under  dynamic  vacuum  for  several  hours  before  usage. 
For  the  synthesis  of  Cu^ySe^Br*,  the  starting  elements  were 
thoroughly  ground,  sealed  in  quartz  ampules,  and  heated  to  1423  K 
for  1  h,  then  cooled  down  to  1073  K  and  annealed  for  48  h.  All 
procedures  were  carried  out  in  box  furnaces  with  heating  and  cooling 
rates  of  5  K/min.  Each  mixture  of  starting  materials  had  the  nominal 
composition  of  Ci^,ySe1_JBrjr  However,  after  the  annealing  step,  a 
small  amount  of  unreacted  elementary  copper  could  be  removed  from 
the  reaction  product,  indicating  that  the  synthesized  samples  are 
slightly  copper  deficient  (Supporting  Information,  Figure  S7).  '  The 
obtained  powders  were  hand  ground  and  consolidated  into  1—1.5  mm 
thick,  12  mm  diameter  discs  at  873  K  for  5  h  under  a  pressure  of  40 
MPa  by  induction  hot  pressing  in  high  density  graphite  dies.19  The 
resulting  pellets  have  more  than  95%  theoretical  density,  determined 
from  the  geometric  densities. 

Characterization.  X-ray  diffraction  experiments  were  performed 
on  a  Siemens  D5000  powder  diffractometer  equipped  with  a  Braun 
M50  position-sensitive  detector,  Ge  (ill)  monochromator  (Huber 
615  002),  and  Cu  radiation,  with  a  step  size  of  0.0078°.  Pawley  fits 
were  performed  with  TOPAS  Academic  V5.0  applying  the 
fundamental  parameter  approach.20  An  exemplary  diffraction  pattern 
including  the  fit  can  be  found  in  Figure  SI.  The  photoelectron  spectra 
(XPS)  were  collected  on  pressed  pellets  with  a  Thermo  Scientific 
ESCALAB  250Xi  with  monochromatic  Al  K„  X-ray  (1486.74  eV)  and 
energy  resolution  of  0.42  eV.  The  surface  was  sputtered  by  a  3  keV 
electron  gun  for  120  s  in  order  to  remove  the  surface  oxide  layer  prior 
to  taking  the  measurement.  The  oxygen  Is  peaks  of  both  samples  were 
completely  removed  after  sputtering  (Figure  S3a).  The  carbon  Is  peak 
(binding  energy  284.5  eV)  was  used  as  a  reference  to  calibrate  the 
binding  energies  of  the  other  core  level  spectra  (Figure  S3b).  The 
relative  contents  of  Cu  to  Se  were  determined  using  a  Thermo- 
Finnigan  ELEMENT2  single  collector  sector  field-inductively  coupled 
plasma  mass  spectrometer  (ICP-MS).  Information  about  sample 
preparation  and  operating  parameters  are  given  in  the  Supporting 
Information.  X-ray  Fluorescence  (XRF)  measurements  were  per¬ 
formed  on  an  18  kW  rotating  anode,  two-circle  diffractometer  with 
monochromated  Mo  radiation  (E  =  17.4  keV)  operating  at  50  kV 
and  160  mA.  Spectra  were  taken  at  an  incident  angle  of  6  =  20°, 
resulting  in  an  average  penetration  depth  of  9.3  ^/m.  The  horizontal  by 
vertical  beam  footprint  at  the  sample  was  0.6  X  6.3  mm,  and  the  total 
collection  time  for  each  sample  was  6  min.  XRF  data  were  collected 
with  a  Vortex-EM  350  //m  thick  silicon  drift  detector.  The  fluorescence 
yield  was  corrected  for  self-absorption  and  detector  efficiency. 
Measured  spectra  can  be  found  in  Figure  S4.  The  thermal  diffusivity 
a  was  measured  with  a  Netzsch  laser  flash  diffusivity  instrument  (LFA 
457)  under  continuous  argon  flow.  In  order  to  maximize  the 


emissivity,  samples  were  spray-coated  with  a  thin  layer  of  graphite 
before  the  measurement.  Thermal  conductivity  was  determined  via  K  = 
aCpd  with  the  heat  capacity  Cp  and  the  geometric  density  d.  Cp  values 
from  previous  publications  were  used,  which  are  in  good  agreement 
with  the  Dulong-Petit  approximation  up  to  600  K,  and  a  Cp 
independent  of  the  Br  content  was  assumed. '21  Electrical  transport 
was  characterized  via  measurements  of  the  Seebeck  coefficient,  Hall 
coefficient,  and  electrical  resistivity.  Electrical  resistivity  and  Hall 
coefficients  were  measured  simultaneously  using  the  van  der  Pauw 
technique  with  pressure-assisted  contacts  under  dynamic  vacuum.  The 
Seebeck  coefficient  was  calculated  from  the  slope  of  the  voltage  versus 
temperature  gradient  measurements  from  chromel-Nb  thermocouples 
under  dynamic  vacuum  with  a  maximum  AT  of  7.5  K  for  all 
temperatures.  Speed-of-sound  measurements  were  performed  at  room 
temperature  to  extract  the  longitudinal  and  transverse  sound  velocities, 
respectively.  The  data  were  obtained  by  use  of  a  Panametrics  NDT 
5800  pulser/receiver  head;  the  response  was  recorded  with  a 
Tektronix  TDS  1012  digital  oscilloscope.  Honey  was  used  as  a 
couplant  between  the  sample  and  the  ultrasonic  transducer. 

■  RESULTS  AND  DISCUSSION 
Chemical  Characterization.  Copper  selenide  has  a  wide 
range  of  stoichiometric  deviation  (Cu2_^Se). 17,22,23  Above  410 
K,  the  compound  has  an  average  cubic  structure  with  Se  ions 
forming  a  face  centered  cubic  lattice  and  copper  ions  being 
distributed  over  tetrahedral  and  trigonal  sites.24  The  room 
temperature  phase  shows  a  very  complex  superstructure.2^ 
Several  structures  were  |>roj3osed  in  the  last  decades  and 
controversially  discussed."  ^  X-ray  powder  diffraction  data  of 
the  synthesized  samples  are  shown  in  Figure  1.  Position  and 
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Figure  1.  Room  temperature  X-ray  powder  diffraction  data  of  the 
synthesized  Cu2.ySe[_^Brx  with  the  formal  composition  Cu2.ySe1_i.BrJ. 
( x  =  0.00,  0.02,  0.03,  0.04,  and  0.05).  A  triclinic  anorthic  unit  cell 
derived  from  LP  search  (TOPAS  Academic  4.1)  was  used  in  order  to 
determine  lattice  parameters  by  Pawley  fits  for  each  composition.  An 
exemplary  Pawley  fit  of  the  Br-free  sample  can  be  found  in  Figure  SI. 
Lattice  parameters  show  no  significant  changes  with  composition; 
nevertheless,  the  obtained  data  are  in  very  good  agreement  with 
literature  diffraction  data  of  Cu2_,,Se.  No  reflections  of  starting 
materials  or  possible  impurity  phases  could  be  observed. 


intensity  of  the  observed  reflections  fit  well  to  other 
experimental  data  of  copper  selenide  in  the  literature. 
Rietveld  refinements  against  standard  laboratory  X-ray  powder 
diffraction  data  with  the  published  structure  models  do  not 
provide  reasonable  information  on  positions  or  site  occupan¬ 
cies.  Se2-  and  Br-  cannot  be  distinguished  by  X-ray  diffraction 
techniques  due  to  their  similar  atomic  scattering  form  factors 
for  X-rays.  However,  lattice  parameters  were  derived  from 
Pawley  fits  using  a  small,  anorthic  unit  cell  derived  from  a  LP 
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Figure  2.  XPS  core  level  spectra  of  Cu2_^Se  (orange)  and  nominal  Cu2_ySe095Br0{15  (blue).  After  sputtering,  only  intensities  from  copper  and 
selenium  can  be  observed  in  the  undoped  sample  and  additionally  bromine  in  the  Br-doped  sample.  Compared  with  standard  CuO  and  Cu20 
spectra,  there  are  no  features  of  oxidation  in  the  investigated  samples.”'  In  combination  with  the  diffraction  data,  the  samples  can  therefore  be 
considered  as  phase  pure.  Comparing  doped  and  undoped  samples  directly,  we  can  observe  slight  changes  in  peak  position  and  peak  shape, 
indicating  that  the  local  bonding  environment  might  be  slightly  different  in  the  investigated  samples. 


search  (TOPAS  Academic  5.0).  Lattice  parameters  for  the 
different  compositions  do  not  vary  significantly  nor  do  these 
values  display  a  trend  upon  increasing  bromide  concentration, 
as  can  be  expected  from  the  similar  ionic  radii  of  Se  and  Br.  The 
standard-laboratory  X-ray  powder  diffraction  data  provide  a 
hint  that  no  secondary  phases  are  present,  as  there  are  no  signs 
of  remaining  starting  materials  or  other  crystalline  secondary 
phases;  however,  it  does  not  fully  corroborate  the  extent  of 
formation  of  a  solid  solution  for  the  substitution  of  Se  with  Br. 
To  check  for  (a)  Br  incorporation,  (b)  possible  CuBr  impurity 
phases,  and/or  (c)  differences  in  the  oxidation  states  resulting 
from  impurity  phases,  X-ray  photoemission  spectroscopy 
(XPS)  was  used  on  pellets  of  the  nominal  compositions 
Cu2_^Se  and  Cu2.ySe095Br005.  Comparing  the  measured  XPS 
spectra  (Figure  2)  of  the  sputtered  Cu2_^Se  sample  with 
previously  reported  XPS  data,34,3S  neither  significant  differences 
in  peak  position  and  peak  shape  nor  additional  peaks  can  be 
observed,  indicating  that  no  side  phases  or  other  oxidation 
states  are  present  in  the  sample.  A  t  higher  values  of  x  (x  > 
0.07),  CuBr  as  secondary  phases  could  be  observed  via  X-ray 
powder  diffraction  (Figure  S2),  and  these  compositions  were 
not  analyzed  in  order  to  ensure  that  no  phase  segregation 
influences  the  transport.  For  copper  selenides,  this  is 
particularly  important  because  they  exhibit  low  carrier 
mobilities,  and  an  impurity  phase  with  high  mobility  would 
dominate  the  transport.''  The  presence  of  bromide  ions  in  the 
doped  sample  can  also  be  verified,  as  the  corresponding  3d 
peak  can  only  be  seen  in  the  doped  sample.  In  addition,  a 
shoulder  on  the  high  binding  energy  side  can  be  observed  in 
the  Cu  2p3/2  and  Se  3d  peaks  in  the  undoped  but  not  in  the 
doped  sample.  The  account  for  the  shoulder  is  not  entirely 
clear,  but  it  might  be  due  to  slight  variations  of  the  local 
bonding  environment  in  the  undoped  sample.  Since  Br-  is 
expected  to  be  an  n-type  dopant  (electron  donor),  additional 
analysis  with  respect  to  the  exact  composition  and  the  Br 
content  is  highly  important  in  order  to  understand  the 
thermoelectric  transport  properties.  X-ray  fluorescence  (XRF) 
analyses  were  utilized  in  order  to  obtain  the  exact  ratios  of 
bromine  to  selenium  in  each  sample.  Because  of  an  additional 
Cu  signal  in  the  XRF  spectra  from  the  sample  holder,  the 
relative  Cu— Se  composition  cannot  be  extracted  via  XRF,  and 
relative  changes  in  the  copper  to  selenium  ratios  were  obtained 
using  inductively  coupled  plasma  mass  spectroscopy  (ICP-MS). 
Assuming  that  the  undoped  sample  has  the  nominal 
composition  of  Cu2Se,  the  changes  in  the  stoichiometry  of 
the  doped  samples  can  be  calculated  and  are  shown  in  Figure  3 


Cu2-ySei_xBrx  -  [Br]  nominal 

Figure  3.  Composition  of  Cll2  ,Se ,  ....Br.  with  the  relative  changes  in 
the  copper  content  [Cu]  determined  using  ICP-MS  and  the  Br 
content  [Br]  determined  by  X-ray  fluorescence  analyses.  With 
increasing  Br  substitution  for  Se,  the  Cu  content  in  the  matrix 
decreases  as  well,  leading  to  a  larger  Cu  nonstoichiometry.  Error  bars 
in  the  ICP-MS  analyses  represent  the  standard  error,  and  the  error 
bars  in  the  XRF  analyses  are  smaller  than  the  symbols. 


and  Table  SI.  First  of  all,  these  data  show  an  increasing  amount 
of  bromide  ions  in  the  samples  with  increasing  synthetic 
content.  While  the  absolute  values  are  slightly  higher  than 
expected,  this  shows  the  stoichiometric  incorporation  of  Br  for 
Se.  Using  the  ICP-MS  data,  the  relative  amount  of  copper  with 
respect  to  selenium  and  bromine  decreases  with  increasing 
bromine  concentration.  In  other  words,  with  the  increasing 
substitution  of  Se  with  xBr  in  Cu^Se^-Br^.,  the  amount  of  Cu 
decreases  as  well,  leading  to  a  larger  Cu  nonstoichiometry  as 
already  observed  for  the  undoped  Cu2_^Se  ( y  >  8). 17,18 
Therefore,  y  is  influenced  by  both  the  intrinsic  Cu  non¬ 
stoichiometry  of  8  in  Cu2_^Se  and  the  change  in  x.  As  the 
absolute  value  of  8  is  small  and  there  is  a  wide  range  of 
stoichiometric  deviation  in  Cu2_^Se,  all  compositions  deter¬ 
mined  for  Cu2.ySe1_vBr1.  are  regarded  with  respect  to  Cu2Se  (<5 
=  0)  in  Figure  3.  Here,  the  relative  changes  in  the  Cu  content 
show  that  when  two  bromide  ions  are  substituted  for  two 
selenium  ions,  the  material  becomes  more  deficient  of  roughly 
one  Cu  ion.  The  occurrence  of  this  interesting  defect  chemistry 
and  its  influence  on  the  thermoelectric  transport  will  be 
discussed  below. 

Potential  for  High  zT  at  Room  Temperature  in 

Cu2_aSe.  Figure  4a  shows  the  temperature  dependence  of 
the  measured  zT  values  of  different  copper  selenide  samples 
with  a  nominal  composition  of  Cu2_^Se.  However,  a  significant 
variation  of  the  measured  transport  can  be  observed,  which  is 
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Figure  4.  (a)  Temperature  dependence  of  the  figure  of  merit  zT  of 
different  copper  selenide  samples  with  a  nominal  composition  of 
Cuj-^Se,  showing  the  spread  of  values  due  to  the  intrinsic  Cu 
deficiency,  (b)  Single  parabolic  band  model  of  the  transport  data 
leading  to  estimates  of  the  effective  mass  m*,  the  mobility  parameter 
//0,  and  the  quality  factor  B  allowing  a  prediction  of  the  thermoelectric 
figure  of  merit  zT  for  different  carrier  concentrations.  The  model 
assumes  that  acoustic  phonon  scattering  limits  the  carrier  mobility. 
Reducing  the  hole  carrier  concentration  in  Cu2_,ySe  by  1  order  of 
magnitude  is  predicted  to  result  in  a  maximum  zT  of  1.16  at  305  K. 


related  to  the  intrinsic  Cu  deficiency  S  in  Cu2_^Se  and  the 
resulting  changes  in  the  Hall  carrier  concentrations.  The  actual 
stoichiometry  of  our  Cu2Se  sample  is  about  Oh 982Se  based  on 
its  transport  properties  in  Figure  5.  A  larger  <5  ultimately  leaves 
holes  h+  as  carriers.  The  individual  components  of  the  figure  of 
merit,  Seebeck  coefficient,  resistivity,  and  thermal  conductivity 
of  one  representative  Cu^Se  sample  can  be  found  in  Figure  5. 
While  the  Cu  vacancy  in  Cu2_^Se  is  of  intrinsic  nature,  the 
holes  behave  extrinsically,  and  a  single  parabolic  band  model 
can  be  applied  to  analyze  the  transport  and  predict  the 
optimum  carrier  densities  for  a  maximum  figure  of  merit.38’39 
The  resulting  prediction  of  zT  versus  Hall  carrier  concentration 
is  shown  in  Figure  4b.  Copper  selenide  is  predicted  to  reach  a 
maximum  zT  of  1.16  at  305  K  with  a  quality  factor  B  of  0.47 
(eq  12).  This  value  is  greater  than  the  state-of-the-art  zT  value 
for  bulk  materials  at  room  temperature,  which  is  around  1.1 
achieved  in  Bi2.xSbvTe3  40'  The  large  zT  value  is  a  result  of  a 
relatively  high  weighed  mobility  (p0m*i/2)  of  around  100 
cm2(Vs)-1,  which  is  comparable  to  the  values  found  for  PbTe 
and  PbSe,"  in  combination  with  a  very  low  lattice  thermal 
conductivity  of  0.16  W(Km)-1.  The  lattice  thermal  conductivity 
was  obtained  using  the  Wiedemann— Franz  law  and  calculated 
Lorenz  numbers,  and  is  much  lower  than  typical  lattice  thermal 
conductivities  of  1  W(Km)-1  found  in  thermoelectric  materials 
and  possibly  a  result  of  the  low  symmetry  crystal  structure  and 
high  Cu  disorder  in  Cu2_^Se  at  room  temperature,  which 
transitions  into  a  molten  Cu  sublattice  at  high  temperatures. 


These  mechanisms  are  also  responsible  for  the  lattice  thermal 
conductivity  of  ~0.22  W(Km)-1  observed  in  Cu7PSe6,  another 
superionic  Cu+  conductor/ 

In  the  case  of  Cu2_gSe,  the  lattice  thermal  conductivity  of 
0.16  W(Km)-1  has  to  be  regarded  as  a  lower  limit  of  K\at  since  it 
has  been  calculated  within  the  single  parabolic  band  model 
under  the  assumption  that  acoustic  phonon  scattering  limits  the 
charge  carrier  mobility.  Differences  in  band  structure  and 
carrier  scattering  assumptions  tend  to  lead  to  an  overestimation 
of  the  Lorenz  number  and  underestimation  of  the  lattice 
thermal  conductivities.4  The  temperature  dependence  of  the 
lattice  thermal  conductivity  below  the  phase  transition  as  well 
as  a  calculation  of  the  glass  limit  for  the  lattice  thermal 
conductivity  can  be  found  in  the  Supporting  Information. 
Modeling  the  thermoelectric  transport  and  predicting  the 
optimum  carrier  density  have  proven  to  be  powerful  tools  in 
thermoelectrics.4'’-46  Figure  4b  shows  that  the  samples  of 
copper  selenide  intrinsically  have  a  Hall  carrier  concentration  of 
4—7  X  102°  carriers  per  cm3,  much  greater  than  the  optimum 
value  of  3  X  1019  cm-3.  Doping  this  material  with  an  electron 
donor  should  reduce  the  amount  of  holes  in  Cu2_^Se  and 
therefore  lower  the  carrier  concentrations  into  a  range  of  much 
higher  figures  of  merit. 

Thermoelectric  Transport  Properties  in  Cu2.ySe,_xBrx. 

Since  copper  selenide  is  intrinsically  a  material  with  extrinsic  p- 
type  behavior,  an  n-type  dopant  (or  electron  donor)  will  be 
needed.  A  possible  candidate  is  Br  as  substituent  on  the  Se  site 
to  form  Cu^Se^-Br^..  Bromine  has  been  shown  to  be  a  very 
effective  n-type  dopant  in  PbTe  and  PbSe  thermoelectrics.47’48 
Copper  selenide  has  been  doped  with  iodine  before,  but  this 
did  not  create  an  appreciable  change  in  the  room-temperature 
zT,  and  no  Hall  carrier  concentrations  were  reported.  ’ 
Measured  transport  properties  of  Cu^ySe^^Brj.  are  shown  in 
Figure  5.  Comparing  the  transport  of  an  undoped  material  ( x  = 
0.0)  with  a  large  hole  carrier  concentration  with  the  data  for 
doped  samples,  it  is  clear  that  substituting  Br  for  Se  has  the 
desired  effect  of  reducing  the  carrier  concentrations.  The 
Seebeck  coefficient  and  Hall  carrier  concentrations  are  all 
positive,  indicating  holes  as  the  majority  carriers  in  this  system. 
The  Hall  carrier  concentration  nH  of  the  holes  decreases  as  the 
amount  of  Br  increases.  The  resistivity  p  and  the  Seebeck 
coefficient  S  both  increase  as  x  increases.  The  samples  with  x 
greater  than  zero  all  have  markedly  lower  thermal  conductivity 
than  Cu2_aSe  due  to  the  reduced  electronic  contribution  via  the 
Wiedemann— Franz  law. 

Defect  Formation  Mechanism  and  Compensating 
Defects.  While  the  substitution  of  Br  for  Se  clearly  decreases 
the  Hall  carrier  concentrations  of  the  holes,  the  actual  changes 
are  lower  than  expected.  The  doping  efficiency  of  Br,  calculated 
with  respect  to  the  high  carrier  concentration  of  the  undoped 
Cu2_^Se  is  below  100%  (Figure  6a).  A  table  of  the  expected  and 
measured  carrier  concentrations  as  well  as  the  doping 
efficiencies  can  be  found  in  Table  S5.  While  Br  usually  shows 
high  doping  efficiencies,30  low  doping  efficiencies  in  thermo¬ 
electric  materials  are  not  uncommon  and  often  occur  visibly  in 
conjunction  with  phase  segregations  due  to  the  Fermi  level 
moving  deeper  into  the  band.  In  this  case,  however,  no 

secondary  phases  can  be  detected  by  powder  X-ray  diffraction 
or  XPS.  Furthermore,  while  usually  the  direction  for  doping  is 
to  move  the  Fermi  level  deeper  into  a  band,  i.e.,  increasing  the 
carrier  concentration,  in  this  case  the  substitution  aims  at 
moving  the  Fermi  level  toward  the  band  edge  in  an  attempt  to 
decrease  the  amount  of  carriers  present.  As  seen  in  Figure  3, 
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Figure  5.  Temperature  dependence  of  the  thermoelectric  transport  properties  of  different  bromine  doped  copper  selenide  samples  Cuj^Sei  ^Br,. 
with  x  =  0.00  (orange  circles),  0.02  (blue  squares),  0.03  (blue  diamonds),  0.04  (blue  upward  triangles),  and  0.05  (blue  downward  triangles). 
Bromine  works  as  an  n-type  dopant  in  copper  selenide.  With  increasing  amounts  of  bromide  ions,  the  electrical  resistivity  p  and  the  Seebeck 
coefficient  S  increase  over  the  investigated  temperature  regime  due  to  the  decreasing  hole  carrier  concentration  nH.  The  total  thermal  conductivity 
decreases  as  well  as  a  result  of  the  reduced  electronic  contribution.  Note  that  the  Hall  mobility  pH  decreases  with  decreasing  Hall  carrier 
concentration,  and  the  zT  values  of  the  doped  samples  only  increase  slightly. 


with  increasing  amounts  of  Br  the  amount  of  Cu  in  the  samples 
decreases.  Starting  from  the  intrinsic  material  Cu2_^Se,  the 
synthesis  does  not  result  in  a  stoichiometric  Cu2Se  but  leads  to 
Cu2_^Se  with  an  excess  amount  of  metallic  copper  after  the 
synthesis  (eq  l).  Even  though  this  excess  amount  of  copper  has 
been  removed  from  the  reaction  product  before  hot-pressing,  it 
can  be  regarded  as  phase  segregation  (eq  l).  The  intrinsic,  Cu- 
defective  Cu2_gSe  can  be  expressed  in  the  Kroger- Vink  notation 
(eq  2), 5  showing  that  the  copper  vacancy  V'Cu  acts  as  an 
acceptor  and  introduces  one  hole  8h‘  per  copper  vacancy.  Note 
that  for  simplicity,  in  eqs  3  and  5  the  reaction  is  written  using 
the  atomic  elements,  not  the  actual  synthetic  reaction  with 
CuBr.  Starting  from  the  intrinsic  material  Cu2_aSe,  the 
synthesis  does  not  result  in  a  stoichiometric  Cu2Se  but  leads 
to  Cu2_^Se  with  an  excess  amount  of  metallic  copper  after  the 
synthesis  (eq  l).  Even  though  this  excess  amount  of  copper  has 
been  removed  from  the  reaction  product  before  hot-pressing,  it 
can  be  regarded  as  phase  segregation  (eq  l).  The  intrinsic,  Cu- 
defective  Cu2_5Se  can  be  expressed  in  the  Kroger- Vink  notation 
(eq  2),51  showing  that  the  copper  vacancy  VCu  acts  as  an 
acceptor  and  introduces  one  hole  8hm  per  copper  vacancy.  Note 
that  for  simplicity,  in  eqs  3  and  5  the  reaction  is  written  using 
the  atomic  elements,  not  the  actual  synthetic  reaction  with 
CuBr. 


Cu2Se  ->  Cu2_sSe  +  <5Cu  (l) 

Cu2_,Se  =  2Cu*u  +  Se*  +  «5V^U  +  «5h*  (2) 

N-type  doping  with  Br-  should  lead  to  a  material  with 
composition  CujSe^Br,.  (eq  4),  if  there  were  no  Cu  vacancy 
formation  involved  at  all,  leading  to  electrons  due  to  the  donor 
Br  on  the  Se  site  Br*Se: 

2Cu  +  (l  —  x)Se  +  xBr  ->  CujSej^Br^.  (3) 

Cu2Sei_ xBrx  =  2Cu£u  +  (l  -  »)Seje  +  xBr*e  +  xe  (4) 

However,  a  Cu  vacancy  mechanism  is  involved  (eqs  5  and  6) 
where  the  resulting  composition  is  Cu^Se^Br,  with  y 
depending  on  the  starting  value  of  8  and  x.  Equation  6)  shows 
in  defect  notation  how  the  electrons,  created  by  Br  donors  Br*Se 
are  partially  compensated  by  the  holes  resulting  from  Cu 
vacancy  formation. 

2Cu  +  (l  —  *)Se  +  xBr  ->  CvL2_ySe1_lfiTx  4-  yCu  (5) 

Cu^Se^Br*  =  2Cu£u  +  yVQn  +  (l  -  *)Sese  +  *Br*e 
+  xe  +  yh*  (6) 

Comparing  the  compositions  in  Figure  3,  it  seems  that  the 
material  becomes  slightly  more  Cu  deficient  with  ~0.5  Cu 
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(1-x)  in  Cu2Sei-xBrx 


Figure  6.  (a)  Doping  efficiency  of  Br  in  Cuj.ySej.^r^.  The  measured 
Hall  carrier  concentration  is  higher  than  the  expected  carrier 
concentration,  calculated  for  one  electron  introduced  per  bromide 
ion  (Table  SS).  (b)  Schematic  defect  formation  enthalpy  for 
Cu2.ySe1_aBrI.  Shifting  the  Fermi  level  to  higher  energies  with 
increasing  Br  content  x  leads  to  a  decreasing  (more  negative) 
formation  enthalpy  of  Cu  vacancies.  The  intrinsic  defect  chemistry 
leads  to  the  observed  lower  doping  efficiency  of  50%,  where  each 
bromide  ion  creates  half  a  killer  vacancy  of  Cu. 


vacancies  per  Br  introduced  into  the  system.  While  bromine 
introduces  an  electron  into  the  system,  the  resulting  0.5  Cu 
vacancies  will  introduce  half  a  hole  in  comparison,  leading  to  a 
doping  efficiency  of  50%.  This  can  indeed  be  observed  for 
lower  values  of  x,  whereas  at  higher  doping  content  the  change 
in  the  carrier  concentration  flattens  out  (Figure  6a).  The  drop 
in  efficiency  can  be  related  to  segregation  (or  nonreaction)  of 
CuBr,  which  should  introduce  even  more  holes.  As  the  samples 
with  larger  synthetic  amount  of  Br  (x  >  0.7)  showed  CuBr  as 
secondary  phases,  this  mechanism  seems  reasonable.  While 
vacancy  formation  has  already  been  shown  to  have  an 
important  effect  on  the  carrier  concentration  in  AZn2Sb2 
Zintl  phases,52’  the  compositional  data  and  the  resulting 
carrier  concentrations  clearly  show  the  influence  of  Cu  vacancy 
formation  as  a  compensating  defect  in  this  material.  Zunger  has 
provided  a  model  to  understand  the  influence  of  defects  on  the 
Fermi  level  in  semiconductors.54  For  instance,  when  an 
electron  producing  dopant  (donor,  Br-  in  this  case)  becomes 
incorporated  in  a  solid,  it  raises  the  Fermi  level.  When  the 
Fermi  level  is  raised  further  into  the  conduction  band,  at  some 
point  the  defect  formation  enthalpy  of  acceptor  defects 
becomes  negative,  and  acceptor  defects  such  as  cation  vacancies 
or  interstitials  are  created.  These  acceptor  defects  act  as  “killer 
defects,”  effectively  pinning  the  Fermi  level  to  a  maximum.5 1  In 
the  case  of  a  p-type  material,  hole  doping  (introducing  an 


acceptor)  will  shift  the  Fermi  level  deeper  into  the  valence 
band,  and  the  competing  mechanism  is  either  the  formation  of 
anion  vacancies  or  cation  interstitials,  which  will  also  result  in 
the  pinning  of  the  Fermi  level  via  killer  defects.  In  the  case  of 
Cu2.ySe!_vBrx,  the  occurrence  of  copper  vacancies  as  “killer 
defects”  occurs  when  moving  the  Fermi  level  closer  to  the 
valence  band  edge  and  not  further  into  the  conduction  band. 
The  intrinsic  material  Cu2_^Se  already  is  deficient  in  Cu,  i.e.,  <5 
>  0,  due  to  entropic  reasons  leading  to  a  temperature  and 
intrinsic  defect  enthalpy  dependence  of  8.  With  increasing  Br 
content,  the  defect  formation  enthalpy  for  copper  vacancies  VCu 
becomes  negative  (Figure  6b,  orange  line).  This  response  can 
also  be  considered  along  the  lines  of  Le  Chatelier’s  principle 
because  adding  Br  alters  a  system  which  has  an  equilibrium  Cu 
defect  concentration  and  Fermi  level.  The  system  then  adjusts 
to  partially  counteract  the  raising  of  the  Fermi  level  caused  by 
bromine  substitutional  defects  by  increasing  the  concentration 
of  the  compensating  copper  vacancy  defect.  In  other  words,  the 
dependence  of  the  defect  formation  enthalpy  for  VCu  in 
Cu2_^Se  is  shifted  to  values  where  any  increase  in  the  Fermi 
level  will  introduce  more  Cu  vacancies,  which  ultimately  act  as 
compensating  defects  in  the  valence  band. 

Charge  Carrier  Scattering  Mechanism.  While  the 
parabolic  band  model  in  Figure  4b  predicts  an  improvement 
of  the  figure  of  merit  at  room  temperature,  the  measured  zT 
values  do  not  exceed  unity,  and  those  at  high  temperatures  are 
only  slightly  greater  than  the  zT  values  of  Cu2_^Se.  This  is  in 
accordance  with  the  findings  for  iodine  doped  Cu2_^Se.5  To 
understand  why  doping  with  Br  did  not  result  in  the  predicted 
high  zT  values,  we  look  at  the  results  of  the  SPB  model  for  data 
at  305  K,  shown  in  Figure  7.  All  parameters  for  the  model  at 
305  K  and  a  full  transport  model  and  description  for  the 
temperature  of  550  K  can  be  found  in  the  Supporting 
Information.  Figure  7b  shows  a  prediction  of  how  S  should 
change  with  nH  in  copper  selenide  using  Pisarenko  relation¬ 
ships.  Orange  symbols  denote  Cu2_aSe,  and  blue  symbols 
denote  Cu2.ySe1_vBra;.  Each  material  was  analyzed  separately. 
The  effective  masses  for  each  material  are  within  2%  of  one 
another,  and  it  appears  that  Br  does  not  change  the  curvature 
and  hence  the  density  of  states  effective  mass  of  the  valence 
band  of  Cu2_aSe.  The  values  for  the  lattice  thermal  conductivity 
(Figure  7c)  are  also  the  same  within  the  measurement 
uncertainty,  0.16  and  0.20  W  (Km)-1  for  Cu2_^Se  and 
Cu^Se^Br,,  respectively.  The  most  striking  difference 
between  the  two  materials  is  the  change  in  the  mobility 
parameter  (Figure  7a).  Copper  selenide  has  a  mobility 
parameter  of  34  cm2(Vs)-1,  but  Cu2.ySe1_i.Bra.  has  a  mobility 
parameter  of  only  17  cm2(Vs)-1.  It  appears  that  adding  Br 
drastically  reduces  the  mobility  of  the  holes  in  the  valence  band. 
In  a  material  influenced  by  a  single  parabolic  band  and  with 
acoustic  phonons  as  the  main  scatterers  of  electrons,  the  Hall 
mobility  should  decrease  with  increasing  Hall  carrier  concen¬ 
trations.55’  The  Hall  mobility  data  for  Cu2_ySe,_xBrx  in  Figure 
7a  exhibit  a  slight  upward  trend  with  Hall  carrier  concentration. 
In  order  to  understand  this  unexpected  behavior,  we  need  to 
incorporate  another  scattering  mechanism  to  understand  this 
contradiction  of  our  assumptions. 

The  scattering  mechanism  for  electrons  and  holes  is 
represented  by  the  scattering  relaxation  time:55’50 


where  e  is  the  dimensionless  charge  carrier  energy.  Each 
scattering  mechanism  has  a  characteristic  equation  for  r0  and 


7023 


DOI:  1 0.1 021  /acs.chemmater.5b02405 
Chem.  Mater.  2015,  27,  7018-7027 


Article 


Chemistry  of  Materials 


8 


n„  /  1020cm'3 


nH  /  1020cm'3 


Figure  7.  Results  for  the  SPB  model  for  Cu^Se  and  Cu2.ySe1_;tBr;c  at  305  K  under  the  assumption  of  acoustic  phonon  scattering.  The  Pisarenko 
relationship  (a)  predicts  how  S  changes  with  nH  in  copper  selenide.  Orange  squares  denote  Cu,_,ySe,  and  blue  circles  denote  Cu2_ySe ,  _;i;Brx.  The 
effective  masses  are  basically  identical,  and  Br  does  not  have  a  huge  influence  on  the  curvature  of  the  valence  band  of  Cu2_,;Se.  However,  the  nH 
dependence  of  the  mobility  parameter  does  not  represent  the  expected  trends  (modeled  lines)  for  acoustic  phonon  scattering  (a)  as  nH  increases 
slightly  with  f.iH  for  Cu^Sej^Br^. 


value  of  A.  When  one  scattering  mechanism  does  not  dominate 
scattering  in  the  material,  multiple  scattering  rates  must  be 
taken  into  account.  Each  scattering  rate  is  equal  to  1/t„  where 
Tj  is  the  relaxation  time  for  a  particular  scattering  mechanism. 
The  scattering  rates  are  treated  like  resistors  in  series,  so  the 
total  scattering  rate  1/r  is  equal  to  the  sum  of  the  individual 
scattering  rates  1/t;  (Matthiessen’s  rule).  The  relaxation  time 
for  acoustic  phonon  scattering  is  given  by1,  ’>s 


JitlAV\d 


-1/2 


V2S  2(wkT)3/2 


(8) 


where  S  is  the  deformation  potential,  v;  is  the  longitudinal 
speed  of  sound,  and  d  is  the  density.  Since  two  kinds  of  atoms 
occupy  the  Se  site  (Se  and  Br),  it  is  logical  to  include  an  alloy 
scattering  mechanism.  The  relaxation  time  for  alloy  scattering 

is59 


8/r 


Tlloy 


—  1/2 


3V2Vatom^(l  -  x)U2; 


\kT ) 


1/2 


(9) 


with  Vatom  =  M{nNAd)~l  the  molecular  weight  is  M,  n  is  the 
number  of  atoms  per  formula  unit,  d  is  the  density,  and  NA  is 
Avogadro’s  constant.  U  is  the  alloy  scattering  potential,  which  is 
similar  to  the  deformation  potential  H  in  that  it  describes  the 
way  the  scattering  mechanism  changes  the  energy  states  of  the 
charge  carriers,  x  is  the  fraction  of  substituted  atoms.  The 
densities  and  respective  speeds  of  sounds  for  the  different 
compositions  can  be  found  in  Table  S2.  As  stated  above,  the 
total  relaxation  time  T  is  related  to  the  total  scattering  rate  1/r: 


-l 

t  = 


1  _i_  -1 

'  Tlloy 


(10) 


The  scattering  mechanism  determines  the  parameter  A, 
which  in  turn  determines  the  order  of  the  Fermi  integrals  used 
to  calculate  the  transport  properties.  Since  Tac  and  TaUoy  have  the 
same  dependence  on  A  (A  =  0),  including  alloy  scattering  into 
the  modeling  does  not  change  the  order  of  the  Fermi  integrals, 
and  can  be  modeled  separately.  Using  r  and  the  effective 
mass  m*,  the  mobility  parameter  fi0  can  be  calculated  via 


—  e-r/m* 


(11) 


Applying  the  combined  acoustic  phonon/alloy  scattering 
model  to  values  of  /u0  computed  for  each  value  of  x  with  the 
deformation  potential  E  and  the  alloy  scattering  potential  U  as 
free  parameters,  we  obtain  a  deformation  potential  of  5.8  eV 
and  an  alloy  scattering  potential  of  1.1  eV.  The  deformation 
potential  obtained  for  Cu2.ySe1_vBrI  is  equal  to  that  obtained 
for  Cu2_^Se,  and  alloy  scattering  potentials  around  1.1  can  be 
found  in  several  other  solid  solutions  like  n-type  PbSe^Te*, 
Cdj.j.Zn^.Te,  or  Al]_xGa,cAs/’0  The  Hall  mobility  versus  Hall 
carrier  concentration  curve  accounting  for  acoustic  phonon  and 
alloy  scattering  is  shown  in  Figure  8,  showing  that  due  to 
increased  alloy  scattering  via  Br  introduction,  a  reduction  of  the 
mobility  can  be  observed.  This  reduction  in  //0  in  the  Br-doped 
samples  due  to  alloy  scattering  reduces  the  quality  factor  B  (eq 
12)  of  Cu2_^Se  from  0.47  to  0.19  (Figure  7d).  The  quality 

factor  B  is  a  measure  of  the  maximum  zT  at  a  given  temperature 
,  .  ,  63,64 

and  is  given  by 
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While  Br  acts  as  an  electron  donor,  i.e.,  it  reduces  the  hole 
carrier  concentration  and  increases  the  Seebeck  coefficient,  the 
additional  scattering  from  the  dissimilar  atoms  on  the  anion  site 
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Figure  8.  Hall  mobility  versus  Hall  carrier  concentration  at  305  K  for 
the  undoped  (orange  squares)  and  doped  (blue)  Cu2.ySe1_aBra. 
samples.  Two  pellets  for  each  composition  with  x  >  0  were  measured. 
The  blue  line  represents  the  model  for  the  Hall  mobility  including 
acoustic  phonon  and  alloy  scattering.  An  error  in  mobility  of  10%  is 
represented  by  the  thick  gray  line.  Alloy  scattering  can  explain 
decreasing  mobility  for  the  doped  samples  using  a  deformation 
potential  of  5.8  eV  and  an  alloy  scattering  potential  of  1.1  eV. 


reduces  the  maximum  possible  zT  to  0.6.  Figure  7d  shows  how 
the  figure  of  merit’s  increase  with  decreasing  carrier 
concentration  slows  down  due  to  a  changing  mobility.  Similar 
effects  have  been  observed  in  lead  chalcogenides  and  has 
recently  systematically  been  investigated  by  Wang  et  al.,50 
where  it  was  shown  that  the  mobility  reduction  due  to  chemical 
substitution  on  a  certain  sublattice  depends  on  its  contribution 
to  the  charge  conducting  band.  The  huge  degree  of  mobility 
reduction  observed  in  Cuj.ySe^Br^.  indicates  a  high  con¬ 
tribution  of  the  anion  atomic  orbitals  to  the  valence  band  states 
close  to  the  Fermi  energy,  which  strongly  influences  the  hole 
mobility.  Because  of  the  complex  crystal  structure  at  room 
temperature,  band  structure  calculations  are  only  available  for 
high  temperature  modification.  However,  the  available 
calculations  show  a  large  contribution  of  the  anion  sublattice 
to  the  valence  band  at  high  temperature.'  It  should  be  noted, 
however,  that  another  possible  scattering  mechanism  could  be 
ionized  impurity  scattering,  which  would  express  itself  as  an 
increase  of  the  Hall  carrier  mobility  at  higher  carrier 
concentrations  as  well,  due  to  screening  effects. 1  ’  One  might 
think  that  ionized  impurity  scattering  would  play  a  role  in 
Cu2.ySe1_vBrx  due  to  the  charge  differences  in  Se2-  and  Br-,  and 
an  upward  trend  of  with  nH  can  indeed  be  seen.  However, 
applying  a  combined  acoustic  phonon  and  ionized  impurity 
scattering  model  to  the  data  at  305  K  indicates  that  ionized 
impurity  scattering  is  not  the  predominant  scattering 
mechanism  in  heavily  doped  Cu2.ySej_j.Brj..  Therefore,  the 
electronic  conduction  can  be  sufficiently  described  by  a 
combined  acoustic  phonon/alloy  scattering  model  for  high 
carrier  concentrations,  which  has  been  successfully  utilized  to 
describe  the  electrical  transport  of  other  heavily  doped  systems 
like  PbTe.20  This  is  encouraging  because  it  suggests  that  doping 
with  something  other  than  Br  may  not  affect  the  mobility 
parameter  possibly  leading  to  larger  quality  factors  and  figures 
of  merit,  whereas  ionized  impurity  scattering  implies  that  any 
dopant  would  reduce  f, i0  and  therefore  reduce  the  maximum 
potential  zT.  In  this  context,  an  investigation  of  corresponding 
cation  doped  copper  selenides  would  be  very  interesting  as  the 
reduction  of  the  mobility  could  be  much  less  pronounced  here 
due  to  a  smaller  contribution  of  the  cation  sublattice  to  the 
valence  band.  To  our  knowledge,  cation  doped  copper  selenide 


has  never  been  successfully  synthesized,  which  might  be  due  to 
the  different  chemistries  of  copper  compared  to  that  of  the 
other  3d  elements.  Our  own  doping  attempts  with  nickel  and 
zinc  led  to  the  formation  of  nickel  and  zinc  selenide  (NiAs 
structure)  even  for  very  low  substitution  grades.  However,  as 
Zunger  has  pointed  out  designing  growth  and  synthesis 
conditions  may  be  effective  in  obtaining  the  predicted  high 
figure  of  merit  without  detrimental  effects  on  the  carrier 
mobility.44 

■  SUMMARY 

We  have  demonstrated  that  Cu2_aSe  has  the  potential  to  be  a 
high-performance  thermoelectric  material  at  room  temperature 
if  the  intrinsic  heavily  doped  character  can  be  altered.  Using 
bromine  as  an  electron  donor,  we  have  successfully  shown  that 
it  is  possible  to  reduce  the  hole  concentration  in  copper 
selenide.  Exploring  the  different  possible  scattering  mechanisms 
for  holes  in  this  material  showed  that  doping  on  the  chemical 
site,  which  mostly  contributes  to  the  valence  band,  leads  to 
enhanced  alloy  scattering  and  a  reduced  mobility  when  doping 
with  bromine.  While  it  shows  that  Br  is  not  an  ideal  dopant  to 
achieve  the  predicted  thermoelectric  performance  in  Cu2_^Se, 
different  strategies  for  further  improvement  should  indeed  be 
possible  because  ionized  impurity  scattering  does  not  play  a 
role  here.  Furthermore,  we  have  shown  the  reduced  doping 
efficiency  of  Br  by  increasing  amounts  of  copper  vacancies  due 
to  Br  substitution  on  the  Se  site.  The  interesting  defect 
chemistry  of  Cu2_^Se  shows,  for  the  first  time,  that 
compensating  defects  and  an  effective  pinning  of  the  Fermi 
level  can  be  achieved  when  shifting  the  chemical  potential 
toward  the  valence  band  edge  as  opposed  to  deeper  into  the 
conduction  band.  This  work  shows  the  need  to  fully  understand 
the  underlying  scattering  mechanisms  and  the  crystal  and  defect 
chemistry  of  a  material  in  order  to  achieve  high  figure  of  merit 
thermoelectrics. 
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